Magnetic and Electric Excitations in Split Ring Resonators 
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We studied the electric and magnetic resonance of U-shaped SRRs. We showed that higher order 
excitation modes exist in both of the electric and magnetic resonances. The nodes in the current 
distribution were found for all the resonance modes. It turns out that the magnetic resonances are 
the modes with odd-number of half-wavelength of the current wave, i.e. A/2, 3A/2 and 5A/2 modes, 
and the electric resonances are modes with integer number of whole- wavelength of current wave, i.e. 
A, 2A and 3A modes. We discussed the electric moment and magnetic moment of the electric and 
magnetic resonances, and their dependence to the length of two parallel side arms. We show that 
the magnetic moment of magnetic resonance vanishes as the length side arms of the SRR reduces 
to zero, i.e. a rod does not give any magnetic moment or magnetic resonance. 

PACS numbers: 160.4760,260.5740 



INTRODUCTION 



The idea of negative index materials (NIMs), i.e. ma- 
terials with both negative electrical permittivity, e, and 
magnetic permeability, /i, was first introduced by Vese- 
lago flf. However, it was only recently that such ma- 
terials were investigated experimentally, 0, 0, 0) S S]- 
Although it has been well known how to obtain e < 
material easily (e.g. using lattice of metallic wires), the 
realization of /i < (especially at high frequencies) re- 
sponse had been a challenge, due to the absence of nat- 
urally occurring magnetic materials with negative /i. In 
1999, Pendry et al. 7] suggested a design made a two 
concentric metallic rings with gaps, called split ring res- 
onators (SRRs), which exhibit a /i < around the mag- 
netic resonance frequency ujm- Immediately after Smith 
et al. [il fabricated the first negative index material 
at GHz frequencies. Recently different groups observed 
0, [13, [m, [l3 indirectly negative /i at THz frequencies. 
In most of the THz experiments, only one layer of SRRs 
were fabricated on a substrate and the transmission, T, 
was measured only for propagation perpendicular to the 
plane of the SRRs, exploiting the coupling of the electric 
field to the magnetic resonance of the SRR via asymme- 
try [l3|. It was realized that one only need the single 
SRR to see the magnetic resonance effects. This way 
is not possible to drive the magnetic permeability neg- 
ative. One reason is that is very difficult to measure 
with the existing topology of SRRs and continuous wires 
both the transmission, T, and reflection, R, along the 
direction parallel to the plane of the SRRs. So there is 
a need for alternative, improved and simplified designs 
that can be easily fabricated and experimentally char- 
acterized. This new design was recently achieved in the 
GHz region [H, [H and the THz region [l6|, |l7l ll8|, lii 
by the use of finite length of wires and the fishnet topol- 



ogy. Very recent work has moved the negative refractive 
index into optical wavelength [23, [2lj . 

In this manuscript we systematically studied the elec- 
tric and the magnetic resonance response of U-shaped 
SRRs for different propagating directions. The effective 
electric permittivity, e, as well as the magnetic permi- 
ability, u will be extracted by the retrieval procedure 
2^ . 23I . 24 1 . In addition, the current distribution along 
the sides of the U-shaped SRR will be numerically cal- 
culated. We show that the magnetic resonance are the 
modes with odd-number of half- wavelength of the current 
density wave, while the electric resonance are modes with 
integer number of whole-wavelength the current density 
wave. In addition we studied the dependence of the elec- 
tric and magnetic resonance as a function of the length of 
the side arms of the U-shaped SRR. It is found that the 
magnetic moment of the U-shaped SRR vanishes as the 
length of side arms of SRR reduces to zero. So there is 
no magnetic moment or magnetic resonance for a metal- 
lic rod. One needs the side arms of the U-shaped SRR 
in order to have a magnetic moment. 



ELECTRIC AND MAGNETIC RESPONSES OF 
SRRS 

A common constituent to provide magnetic response 
in metamaterials is the Split-Ring resonator (SRR). The 
SRR in its simplest form consists of a highly conductive 
metallic ring which is broken in one (or several) loca- 
tion(s) by a non-conductive gap of air or other dielectric 
materials. If this ring is placed in a temporally varying 
magnetic field an electric circular current is induced in 
the metallic ring which in turn leads to charge accumu- 
lating across the gaps. The electric field which builds due 
to the charge at the gap counteracts the circular current 
leading to energy stored (predominantly) in vicinity of 



2 



the gaps and magnetic field energy concentrated in the 
region enclosed by the ring. The SRR is thus a resonator 
which couples to a perpendicular magnetic field and can 
be characterized by the effective capacitance of the gaps 
and effective inductance of the loop define by the ring. It 
can be understood in terms of a resonant LC circuit with 
a resonance frequency = 1 /iC, where L is the induc- 
tance and C is the capacitance of the SRR. The resonant 
response of the circular current in the SRR to an exter- 
nal magnetic field leads to a resonant magnetic moment 
which may reach large negative values for array of SRRs 
such that the size of the SRR is much smaller than the 
wavelength of an incident electromagnetic wave around 
the resonance frequency behaves as a homogeneous ef- 
fective medium with at negative (resonant) permittivity 
Meff(w). 

In this letter we report numerical results of electric 
and magnetic responses of a periodic lattice of SRRs for 
different orientations of the SRRs with respect to the ex- 
ternal electric field, -E, and the direction of propagation, 
k. 





FIG. 1: (a) Schematic of periodically arranged U-shaped SRR 
arrays; (b) a single unit cell with geometrical parameters; 
(c) three different configurations of incident electric field, E, 
magnetic field, H , and wave vector, k. 

Our numerical simulations were done with CST 
Microwave Studio (Computer Simulation Technology 
GmbH, Darmstadt, Germany), which uses a finite- 
integration technique, and Comsol Multiphysics, which 
uses a frequency domain finite element method. The 
schematic of the periodic U-shaped SRRs arrays and the 
geometry of a single unit cell used in our numerical simu- 
lation were shown in Fig. [Ija) and (b), respectively. The 
SRRs exhibit different responses to the incident electro- 
magnetic (EM) wave with respect to different configu- 
rations of incident electric field, magnetic field, H, 
and wave vector, k, as shown in Fig. [ijc). First, it is 
well know that the incident EM wave excites a magnetic 
resonance at ujm [1, S 0, S S 01 if the external mag- 
netic field, H , is perpendicular to the SRR plane and the 
wave vector, fc, is parallel to the SRR plane (Fig. [TJc.i)). 
Second, farther study reveals that the incident EM wave 
with k perpendicular to the SRR plane and E parallel 
to the bottom part of SRR (Fig. [Hc.ii)) also excites a 
magnetic resonance at Wm ■ This magnetic resonance 



results from a circular current induced by the external 
electric field, E, because of the asymmetry of the SRR in 
the direction of E, and therefore was called the electric 
excitation coupling to the magnetic resonance (EEMR). 
The EEMR is very valuable for the experimental demon- 
stration of the magnetic resonance at optic frequencies, 
because it is very difficult to measure the transmission 
and reflection with the incident EM wave parallel to the 
SRR plane at high frequencies [i, 0, [Til, [ll IH] . Third, 
the SRRs also exhibit a short-wire-like electric resonance 
[2^ at uo with the incident electric field, E, parallel to 
the side part of SRRs (Fig. [Ijc.iii)). The electric res- 
onance frequency ujq depends on the length of the side 
part of SRRs, ly, being higher for shorter ly. 
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FIG. 2: Extracted effective permittivity Re(e((.<j)) (a) and ef- 
fective permeability Re(^(cij)) (b) for the magnetic resonance 
(red solid), the electric excitation coupling to the magnetic 
resonance (EEMR) (blue dashed) and the short-wire-like res- 
onance (green dotted) of USRRs. The geometric parameters 



are ~ a.y ~ ^ /im, — 
= 0.1 /im and t — 0.05 fim. 



0.2 fim, Ix ~ ly 



0.8 fim, 



24 |. we calcu- 



Using the retrieval procedure [22l . \2c 
lated the effective permittivity e(a;) and permeability 
/i(a;), both real and imaginary part, from the simulated 
transmission, T, and reflection, R. Fig. [2] shows the ex- 
tracted real part of the effective permittivity Re(e([j)) 
and permeability Re(^(tt')) for three different polarized 
incident EM waves as shown in Fig. [TJc). As expected, 
the magnetic resonance, measured by /i (red solid)), oc- 
curs at around uJm — 55 THz. The EEMR gives a sim- 
ilar resonance at roughly the same frequency, shown as 
a resonance shape in e (blue dashed) , which indicate the 
response is due to the electric coupling of the incident 



EM wave to SRR. At a higher frequency, ujq = 135 THz, 
the short-wire-hke electric resonance of the SRR occurs, 
given by the resonance behavior of e (green dotted). Due 
to the periodicity effect, whenever a magnetic resonance 
occurs in Re(/i), an electric anti-resonance will be seen in 
Re(e) simultaneously and vice versa 24i Hi- The period- 



icity effect becomes more noticeable as the wavelength of 
EM wave is comparable with the lattice constant of the 
SRR lattice. For the magnetic resonance (red solid), the 
lattice constant, — i Jim, in the propagating direction 
along the wave vector, fc, is comparable to the resonance 
wavelength. Am = 5.45 /^m, hence, a strong distortion 
in the negative part of resonance in Re(/i) and a signifi- 
cant anti-resonance in Re(e) were observed. On the other 
hand, for the EEMR, the lattice constant, Ux = 0.2 ^m, is 
much smaller than A™, so the periodicity effect is much 
weaker and therefore a sharp resonance in Re(e) (blue 
dashed) with a weak anti-resonance in Re(/i) occurs. 
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HIGHER ORDER EXCITATION MODES OF 
SRRS 

Beside the typical electric and magnetic response, the 
SRRs also exhibits higher order excitation modes. As 
shown in Fig. [3lja), for the incident EM wave with 
the propagating direction perpendicular to the plane of 
SRRs, we found three successive dips at 55 THz, 155 THz 
and 300 THz in the transmission spectra. Correspond- 
ingly, at the same frequencies, we found three resonance 
shapes in Re(e) (Fig. [3jc)). The first excitation mode 
occurs at the lowest frequency fmo = 55 THz, which is 
the usual EEMR response, and has the strongest reso- 
nance in Re(e). The magnitudes of resonance of the sec- 
ond and third excitation modes, occur approximately at 
3/mO and B/mO respectively, and are much weaker then 
the first excitation mode. Higher order excitation modes 
of the short-wire-like response also been observed (Fig. 
[3Kb) and (d)). In Fig. [3] (b), four dips were found in 
the transmission spectra, which result from the first, sec- 
ond and third order electric excitation modes of SRRs. 
Detailed studies show that the second (205 THz) and 
the third (265 THz) modes have similar current density 
distribution, so they are both considered as the second 
excitation mode. 

As a very crude picture we could imagine the reso- 
nances of the SRRs as charge density waves on a rod of a 
length equal to the arc length of the SRR ring. This rod 
supports plasmonic modes [27| (plasmonic means in this 
situation that the inductance is coming from the electron 
mass and the capacitance from the external electric field 
over the surface of the rod) which occur at discrete fre- 
quencies whenever we have current nodes at the ends of 
the rod. This picture is of limited value for the following 
reasons: (i) the EM response of the modes, especially the 
classification of electric vs. magnetic response, depends 



FIG. 3: Transmission spectra (a,b) and the extracted per- 
mittivity, Re(e), (c,d) of the U-shaped SRRs response to the 
incident EM wave. The directions oi E, H and k were shown 
in the insets. 



on the geometry and is entirely different for a U-shaped 
SRR than for a straight rod although the qualitative cur- 
rent distribution over the arc length is equivalent; (ii) the 
plasmonic dispersion depends on the geometric induc- 
tance, i.e. the energy stored in the external field outside 
the metal, which is much larger for the SRR than for 
the rod and spatially non-uniform along the SRR ring; 
and (iii) the modes of coupling are different. Neverthe- 
less, considering the current distribution is an essential 
tool for understanding the resonant modes of a SRR at 
higher THz frequencies. In the microwave region only 
the remnants of the lowest order mode survive. 

Fig. IHshows the distribution of the current and charge 
density for the lowest three resonant modes of a U-shaped 
SRR of side length, ly — 800 nm. The current density is 
obtained directly from the simulation results, the charge 
density relates to the perpendicular electric field at the 
surface of the metal. All fields are time harmonic; the 
current distributions are shown temporally 7r/2 phase 
shifted against the charge distribution. The first three 
panels (a,b,c) of Fig. [4] show the lowest three EEMR res- 
onances for normal incidence to the SRR with the elec- 
tric field breaking the symmetry of the SRR and thus 
we have coupling to the "magnetic" resonance [l3(. The 
surface electric field distribution is qualitatively equiv- 
alent to the pure magnetic coupling (i.e. propagation in 
the SRR plane with perpendicular magnetic field), which 
case is however hard to realize experimentally at such 
high THz frequencies. All three modes have non-zero 
magnetic moment coming from all three "arms" of the 
SRR. The number of current nodes (where charge ac- 
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FIG. 4: Distribution of the perpendicular component of the 
surface electric field (color scale; red positive, blue negative) 
and the bulk current density (arrows) for the lowest few res- 
onant modes of the SRR. The SRR metal is made of a Drude 
model for gold {fp = 2175 THz, = 6.5 THz), the geometry 
parameters are: Ux = ay = 1 nm, — 200 nm (unit cell 
size) , Ix = ly = 800 nm (arm length) , w = 100 nm, t = 50 nm 
(ring width and thickness, respectively). The current distri- 
butions are shown temporally 7r/2 phase shifted against the 
charge distribution. 



cumulates) increases with the resonance frequency from 
two (fundamental mode, corresponding to the resonance 
of the effective LC circuit), one on either side of the 
"gap", to four and six, which have additional nodes in- 
side the continuous metal. The electric excitation of these 
modes occurs via the polarization of the bottom arm of 
the SRR by the electric field of the incident EM wave. 
For a straight rod, these SRR modes would correspond to 
the A/2, 3A/2, and 5A/2 mode (Fig. E^a)); however, for 
the metallic rod, there is no magnetic moment associated 
with the excitation modes. Note that these modes also 
possess an electric dipole moment; therefore the SRR has 
a combined magnetic and electric response in this con- 
figuration. For the other polarization, normal incidence 
to the SRR with the electric field along the symmetry 
axis of the SRR, shown in panels (d,e,f) of Fig. [H we see 
the analog plasmonic modes with three, five and seven 
current nodes; all of which reflect the mirror symmetry 
of the SRR and can thus have no magnetic moment. For 
the metallic rod, they would correspond to the even. A, 
2A, and 3 A modes (which cannot be exited for the rod 
because of their vanishing dipole moment). In the case 
of the SRR they do possess electric dipole moment in the 
E-direction, i.e. along the two parallel arms of the SRR, 
and represent a purely electric response of the SRR. 

Fig. ^a) and (b) shows the distribution of the to- 
tal current (current density j integrated over the cross- 
section of the SRR ring) over the arc length around the 
SRR ring for the three resonant modes of Fig. ma,b,c) 
and Fig. [ll^d,e,f ) respectively. We clearly see the differ- 
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FIG. 5: Current distribution of the lowest 3 modes, (a) A/2, 
3A/2, and 5A/2 mode for the EEMR response; (b) A/2, 3A/2, 
and 5 A/2 mode for the short-wire-like response. Due to the 
nonzero response of lower order modes, nodes of higher order 
modes only reach zero at the positions of the lowest nodes. 



ent number of current nodes, which do not reach zero (for 
the higher modes) because of the superposition with the 
non-resonant response of the lower order modes (which 
are fairly broad due to the high losses). Also note the 
non-uniform spacing of the nodes for the higher orders 
which are different from the straight rod and are due to 
the curvature of the SRR and the coupling to the bottom 
arm or side arms only. 



ELECTRIC AND MAGNETIC MOMENTS OF 
SRRS 

We also studied the electric moment and the magnetic 
moment for the electric excitation of the magnetic reso- 
nance (EEMR) and the short-wire-like resonance of the 
U-shaped SRRs with different length of two parallel side 
arms, ly. The electric moment is calculated by the inte- 
gration of the subtraction between the electric displace- 
ment D and vacuum electric displacement eoE over the 
volume of the whole unit cells. 



J {D-eoE)dr 



Then the electric polarization, normalized by incident 
electric field, Eq, P = p/ (VEq) is shown in Fig. [Sfa) and 
(b) for the EEMR and the short-wire-like response of the 
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FIG. 6: Magnitude of the normalized polarization density, 
P, of the U-shaped SRRs with ly = 0.8 /im (red solid), 0.4 
fj,m (blue dashed) and 0.11 /im (green dotted), respectively, 
(a) Px component as E parallel to the bottom part of SRRs. 
The other two components Py and Pz are nearly zero; (b) Py 
component as E parallel to the side part of SRRs. The other 
two components P^ and P^ are nearly zero. The polarizations 
of the incident EM wave are shown as the insets in the panel 
(a) and (b). 



SRRs with different lengths of the two parallel arms, ly. 
It is clearly seen that the lowest order magnetic reso- 
nance provides the strongest magnetic response; but also 
the higher modes have non-zero magnetic moment. As 
one expects, reducing the length of the parallel arms re- 
duces the magnetic moment. It also shifts the magnetic 
resonance to higher frequencies. This is a combined ef- 
fect of shortened arc length and reduced the geometric 
loop inductance. In the limit of only the bottom arm left 
(i.e. a straight rod) the magnetic response would vanish 
(green dotted curve), and no magnetic response is seen, 
as expected. In this limiting case, the resonance in Re(e) 
still exists, which is a result from the electric resonance 
of the bottom arm of SRR. Note that the magnitude of 
the resonance in Re(e) decreases much slower than the 
magnetic moment does as the arm length, ly^ decreases. 
The reason is that the magnetic moment, m^, resulting 
from the circular current flowing in the loops of the SRR, 
vanishes immediately as the side arms disappear, but the 
resonance in Re(e) always exist as long as the short-wire- 
like electric resonance exists in the bottom part of the 
SRR. However, as shown in green curve of Fig. Elja), 
this short- wire- like electric resonance can not provide any 
magnetic moment. 

The short-wire-like resonance with the propagating di- 
rection of the incident EM wave perpendicular to the 
SRR plane and the electirc field, E, parallel to the side 
arms of the SRR (Fig. [TJcii)), has zero magnetic mo- 
ment. 



SRRs respectively. It is clearly seen that the lowest mode 
of resonance has the strongest electric moment although 
the higher order excitation also has nonzero value. In Fig. 
inila), the electric polarization, Pj, decreases as the length 
of two parallel arms, ly, decreases, because the EEMR 
response becomes weaker and finally disappeared as the 
length ly close to zero (green dotted curve). However, the 
electric polarization, P^, does not vanish in the limiting 
case of Zj, = 0, (green dotted curve), because the EEMR 
degenerates to the short-wire-like electric resonance of 
the bottom arm of the SRR which results in the non- 
zero electric polarization. The polarization density along 
the y-dircction, Py, of the short-wire-like resonance of the 
SRR (Fig. m^b)) with the incident electric field E parallel 
to the side arms of the SRR, also decreases and shifts to 
higher frequencies as ly decreases. 

The magnetic moment is calculated by employing the 
formula: 

J fx (D- eoE)df 

In Fig. [71^ a), we show the magnetic moment, normal- 
ized by the magnetic field of the incident EM wave, Hq, 
— rnz/{VHo), and in Fig.jT^b), the extracted permit- 
tivity, Re(e), for the magnetic resonant modes of the SRR 
as a function of frequency for three different U-shaped 



CONCLUSIONS 

We systematically studied the electric and magnetic 
resonances of U-shaped SRRs with respective to different 
polarizations of the incident EM wave. Higher order ex- 
citation modes were found in both electric and magnetic 
resonances. We show that the magnetic resonance are the 
modes with odd- number of half- wavelength of the current 
density wave, while the electric resonance are modes with 
integer number of whole-wavelength of the current den- 
sity wave. In addition, the current density distribution 
of the lowest three excitation modes were given. We also 
studied the magnetization density as a function of the 
length of the side arms of the U-shaped SRRs. It turns 
out that the magnetic resonance vanishes as the length 
of side arms reduces to zero, i.e. a rod does not give any 
magnetic moment or magnetic resonance. 
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FIG. 7: (a) Magnitude of the normalized magnetization, Mz, [16] 
and (b) the extracted permittivity, Re(e), of the U-shaped 
SRRs with the length ly = 0.8 /xm (red solid), 0.4 jj.m (blue 
dashed) and 0.11 fim (green dotted), respectively. The po- 
larizations of the incident EM wave are shown as the insets [17] 
in the panel (a). The short-wire-like resonance with incident 
EM wave polarized as shown in Fig. [l][c.ii) has zero magnetic 
moment, and therefore is not shown here. 
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